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We show that for a generic quantum mechanical system with more than one open scattering 
channel, it is not possible to fully reconstruct the theory's S'-matrix from spectral information 
obtained in large finite volumes with periodic boundary conditions. Physically distinct S-matrices 
can have identical finite-volume spectra for large finite boxes of arbitrary sizes. If the theory is 
not time-reversal symmetric, there exists an uncountably infinite set of distinct ^-matrices with the 
same spectra. If the theory respects time-reversal symmetry there exists a discrete set of ^-matrices 
with identical energy levels for finite boxes. We illustrate the issue for simple quantum mechanical 
systems in 1+1 dimensions. 



I. INTRODUCTION 

Because it is strongly coupled at low energy, quantum 
chromodynamics (QCD) is not amenable to perturbative 
techniques there. Fortunately, lattice simulations pro- 
vide a first-principles technique for extracting predictions 
from the theory. The lattice has achieved unprecedented 
accuracy for many QCD observables, including some that 
arc difficult to extract experimentally [1]. Indeed, it has 
become a precision tool for many purposes. However, 
not all QCD observables are easily accessible via lattice 
simulations. One important class of these are scattering 
observables in regimes where there is more than one open 
channel. 

Before early groundbreaking work [2-7] it was not gen- 
erally understood that at least some types of scattering 
information are accessible via finite- volume lattice stud- 
ies, which by construction always yield a discrete spec- 
trum of energy levels. Now it is well-understood that in 
cases where the kinematics restrict the scattering to a 
single two-body scattering channel, one can extract the 
phase-shift for the channel given knowledge of the energy 
levels in a finite box. In particular given a sufficiently 
large box, knowledge of the size of the box, and the en- 
ergy eigenlevels in that box, one can extract phase-shifts 
for those energies. By scanning through different box 
sizes — and thus different energies — one can extract the 
phase shift as a function of energy. This approach has 
already proven to be practical for cases of low-energy 
scattering [8-11]. However, this technique, usually re- 
ferred to as "Liischer's formula" , only provides access 
to certain kinds of data. It does not apply in its origi- 
nal form to situations in which there is more than one 
accessible channel. 

This difficulty is acute for the problem of hadron spec- 
troscopy, a subject in which there has been considerable 
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interest from the lattice community [9, 12-22]. The diffi- 
culty is that the vast majority of hadrons are resonances 
above the three-body threshold for which the original 
formulation of the Luscher formula is inapplicable. The 
current state of the art for high-lying resonances is sim- 
ply to calculate in a given relatively small box size and to 
identify the resonance's mass with an energy level. Such 
a procedure is presumably a valid way to identify the 
position of the resonance in the limit of an infinitely nar- 
row resonance and in a regime where the box size is large 
enough to avoid large finite volume effects. However, the 
hadrons of physical interest typically have fairly substan- 
tial widths so this procedure is questionable. Moreover, 
the "mass" of a resonance is a quantity which depends on 
how one separates peaks from background in scattering 
data; such separation inherently creates model depen- 
dency. It would be far preferable to have access directly 
to the scattering observables in the form of the S'-matrix 
directly from the lattice. This would enable one to make 
contact with the physics associated with hadronic widths 
and branching ratios. 

The need for going beyond a simple association of a 
hadron's mass with the energy of a resonance is particu- 
larly clear in cases where the width of a resonance is large 
and multiple channels are involved. A good example is 
the Roper resonance: although it has a low excitation en- 
ergy, it is quite broad and is strongly coupled to the nu- 
cleon plus two pion channel. The Roper resonance is sur- 
prisingly difficult to characterize experimentally and the 
width and branching ratios are poorly determined [23]. 
Thus for the Roper resonance, where an ab initio calcu- 
lation of the scattering observables would be particularly 
helpful in trying to understand the nature of the state, a 
direct extraction of the mass from the energy on a single 
finite volume lattice is particularly difficult. 

There are other cases where a direct extraction of the 
S'-matrix is important. For example, in the on-going pro- 
gram to attempt to understand nuclear phenomena from 
QCD, low energy nuclear scattering observables are of 
interest. The 1 So channel in nucleon-nucleon scattering 
involves a single channel; it's S- matrix is fully character- 
ized by a phase shift and the Luscher formula directly 
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applies. However, scattering in the deuteron channel is 
coupled and includes mixing between the 3 Si and 3 Di 
partial waves and the Liischer formula does not directly 
apply. To proceed with the program of extracting nuclear 
physics from lattice QCD it is important to understand 
how to extract the S-matrix from the lattice. 

There have been significant recent theoretical develop- 
ments aimed at generalizing the Liischer formula relating 
the energy levels in a finite box with periodic boundary 
conditions to many-channel S- matrices [24-27]. Underly- 
ing these works is the key insight behind the Liischer for- 
mula: when the box is significantly larger than the range 
of interactions between the particles, then the particles 
are essentially non-interacting when they encounter the 
boundary conditions associated with the box. Thus the 
S'-matrix alone determines whether the boundary condi- 
tions are met. One expects that complete knowledge of 
the S'-matrix will determine the energy levels in the box. 

This paper concentrates on the inverse problem: If one 
has access to the energy levels of a system computed on a 
set of finite boxes with periodic boundary conditions, can 
one reconstruct the full multi-channel S-matrix? This 
is the nature of the problem presented by lattice QCD. 
In the case of Liischer's formula, in which only a single 
channel is relevant, this is straightforward. The existence 
of a bound state at a given energy is fixed by a single 
piece of information — the phase shift at that energy — 
which can be read off directly given knowledge of the 
box size. By working at a variety of box sizes which 
act to change the energies, one can then map out the 
phase shift as a function of energy. This is tantamount 
to learning the entire S-matrix. However in the multi- 
channel case, the S-matrix contains more than a single 
piece of information at one energy. It is thus not possible 
to determine a multi-channel S-matrix at a given energy 
from the knowledge of levels in a single finite box — even 
if the discrete energy levels in the box include the energy 
of interest. 

Given this situation, it is immediately apparent that 
if it is possible to reconstruct the full S-matrix at all, it 
will necessarily involve extracting the energy levels from 
more than a single box. However, it is not immediately 
obvious whether it is possible to extract the S-matrix at 
a given energy even from the knowledge of energy levels 
in many different boxes. Note that, in the case of a single 
channel problem where the Liischer formula holds, if one 
has two different box sizes which share a common energy 
level, they give redundant information: both boxes allow 
for the extraction of the same phase shift at that energy 
(up to finite volume corrections) . It is by no means clear 
that in a multi-channel situation, where many parame- 
ters arc needed specify the S-matrix, that different box 
sizes will also not give some redundant information in 
such a way as to prevent one from fully recovering the 
S-matrix. Indeed we will show that in general it is not 
possible, even in principle, to fully reconstruct the /S'- 
matrix at some given energy, even with a knowledge of 
the spectrum for an arbitrarily large number of distinct 



box sizes. In the general case, there exist physically dis- 
tinct S-matrices which yield the same same spectra in 
any finite box subject to periodic boundary conditions. 

This paper is aimed at clarifying issues of principle. 
Given this, it is useful to study the simplest context for 
which the issues arise. To do this we will make several im- 
portant departures from the case of greatest interest — a 
system with 3+1 space-time dimensions, with fully rela- 
tivistic kinematics and with any number of particles in 
the final state. 

Instead, we confine our detailed studies to systems with 
I+I space-time dimensions. The full 3+1 dimensional 
problem is much more complicated technically. This 
complication arises from the fact that the free space the- 
ory is rotationally invariant, giving rise to a partial wave 
decomposition for the S-matrix, while the finite box only 
has a discrete subgroup of the rotations associated with 
the symmetries of a cube. However, it seems quite likely 
that the answers to the most fundamental questions of 
principle are not likely to depend on the number of space- 
time dimensions. It seems almost certain that if it turns 
out not to be possible to extract S-matrices from the 
finite volume data for 1+1 dimensional systems then it 
will not be possible for 3+1 dimensional systems. Thus, 
as a first step it is seems sensible to illustrate the im- 
possibility of specifying the S-matrix using a simple 1+1 
dimensional system. 

In considering multi-channel S-matrices there are two 
qualitatively different kinematic regimes to consider. In 
one case, all of the channels involve two particles in the 
final state. In the case of QCD one might consider scat- 
tering between a pion and a strange baryon (A or E) in 
the I — 1 channel in the regime where both the 7r-A and 
7T-S channels are kinematically accessible but below the 
threshold for the ir-ir-A channel. Thus, while there are 
multiple channels in such cases the S-matrix at any fixed 
energy can be fully specified by a discrete set of parame- 
ters. The other case is when the system is in a kinemat- 
ical regime where there may be more than two particles 
in the final state. In this second case one needs to spec- 
ify continuous functions to fully describe the S-matrix 
at a fixed energy. This second case is clearly more chal- 
lenging. Again, our aim here is to keep things as simple 
as possible to illuminate the underlying issues. Thus, we 
will focus in this work on the first case — multiple channel 
problems with only two-particle final states. 

Finally, the analysis of any particular problem depends 
on kinematical details. These simplify for the case of 
nonrelativistic kinematics. In the models studied in de- 
tail here we will restrict our attention to nonrelativistic 
systems, but our result easily generalizes to relativistic 
circumstances. 
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II. A GENERALIZED LUSCHER FORMULA 
FOR A SIMPLE TOY PROBLEM 

Consider the following nonrelativistic two-channel 
scattering problem in 1+1 space-time dimensions. Sup- 
pose that there are two possible channels; each one con- 
sists of two identical bosons of mass M but the parti- 
cles in each of the two channels are different. These bo- 
son should be considered as composite objects with finite 
size — as mesons are in QCD. The energy in channel A is 
simply the kinetic energy while the energy in channel B 
is offset by an amount A. One can view A as the amount 
of energy needed to excite the two particles of type of A 
into those of type B. Since for simplicity the dynamics 
are nonrelativistic in both channels the excitation energy 
does not effect the mass: it must be conserved and thus 
is the same in both channels. As will be discussed later, 
the restriction to non-relativistic dynamics is inessential. 
The scattering problem consists of firing two particles of 
either type A or type B at each other and observing the 
outgoing particles which are either of the original type 
or through the interaction both converted to the other 
type. This problem does not correspond to any interest- 
ing physical system. However, it is a useful toy model for 
illustrating the relevant physics. 

First consider the kinematic restrictions on S. By 
Galilean invariance we expect the 5-matrix that de- 
scribes this physical process to conserve energy and mo- 
mentum, and moreover to be independent of the total 
momentum. We denote momentum of the two particles 
of species j = A,B by kj-\$ where 1,2 labels particles 1 
and 2 in the channel. The total momentum K is obvi- 
ously 



where Kj = — kj-p) is the relative momentum of the 

two particles of species j, and fi = M/2 is the reduced 
mass. The asymptotic states that S is kincmatically al- 
lowed to connect, therefore, have relative momenta re- 
lated by 



4 = k 2 - 2/jA 



(4) 



where we use the shorthand Ajk = Aj — A&. 

Let us assume that the box is much larger than the 
range of the interaction between the particles. Then, in 
large spatial regions the energy eigenstates in the box 
will approximate well the asymptotic form. Because we 
are interested in periodic boundary conditions, outgoing 
asymptotic states will reach the box boundary and be- 
come incoming states on the opposite side of the box. 
The energy eigenstatc are asymptotically of the form 



> asy = E 



j,K jt K) (5) 



3=A,B 



+dj | out; j, Kj, K ) 



where the related to E by Eq. (3) and "in" and 

"out" label the incoming and outgoing amplitudes re- 
spectively. 

Now let us consider the action of S on asymptotic 
states. It relates the coefficients by 



K = kj-i + kj- 



r,2- 



(1) 



Since we are considering a nonrelativistic problem, the 
energy can be written 



1 

2M 



E= — {kl 1 +kl 2 ] 



Aj (2) 



where are Aj = for j = A and Aj = A for j = B. It is 
easy to see that 



E 



K 2 



2(2M) 2fi 



(3) 



S(E,K) jk a k . 



(6) 



Note the S depends only on the relative momenta; the 
dependence on E and K should be understood to be en- 
tirely through Eq. (3). 

Of course, only wavefunctions with certain energies will 
fit in a box of given size, and these energies will be de- 
termined by the dynamics encapsulated in S. In posi- 
tion space, the pieces of this asymptotic wavefunction 
for channel j are 



j, X!,x 2 | in j, Kj, K 



j, X\,X 2 \ OUt j, Kj, K 



asy 



e iKX 

AKX 



asy v^Tm 
I — 



6{x)) 

(e iK J x 9{x) + e- lK i x 9(-x)) 



(7) 



where X is the location of the center of mass \{x\ + x%) and x = (X1—X2) is the relative coordinate, 9 is the usual 
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Heaviside step function, and the normalization is such 
that the incoming probability flux and outgoing proba- 
bility flux are the same if S is unitary. Here, "asymp- 
totic" is understood to mean that the magnitude of the 
relative coordinate x is much larger than range of the in- 
teraction. The symmetric form reflects the fact that the 
particles are identical bosons. 

With these explicit asymptotic wavefunctions in hand, 
we note that Eq. (3) does not fix the sign of Kj. Either 
sign is acceptable since it can be reabsorbed in the def- 
initions of aj and ay; by convention we define Kj to be 
positive. 

Typical lattice simulations are performed with either 
periodic or twisted boundary conditions. So far as the 
issues of principle addressed in this paper are concerned, 
one may verify that twisted conditions offer no additional 
leverage for the problem of extracting S in this toy prob- 
lem. In any event, they may be included with minimal 
effort. For simplicity we restrict our attention to simple 
periodic boundary conditions. If the box under consider- 
ation has spatial extent L, periodic boundary conditions 
imply that the wavefunction obeys 



L 



■r-2 



9(E,K) ) = 



L 



■r-2 



(8) 



In particular, we can project this equality onto each chan- 
nel and require the equality to hold separately for the 
pieces which depend on exp(ifc j; iX2) and on exp^/c^a^)- 
The kj-i -dependent equality implies 



while the fcj^-dependent equality implies 



a : = a j e 



-ikj-iL 



(9) 



(10) 



By comparing these two equations we see that all chan- 
nels imply the same quantization condition, 



KL = 2itZ ; 



(11) 



the integer Z characterizes the total momentum. 

Note that in lattice simulations one can restrict at- 
tention to correlation functions of operators with a fixed 
known momentum and thus extract energy levels with a 
fixed known Z. We will assume that this can be done 
for this toy problem as well. In principle one can try 
extracting the S'-matrix using a set of measurements at 
any fixed Z or a mixed set of Zs. Again, as far as issues 
of principle are concerned, this freedom gives no addi- 
tional power for extracting information about S. We will 
henceforth assume Z is fixed and suppress dependence 
on the center of mass momentum K. 

The relations between the coefficients of asymptotic 
incoming and outgoing wavefunction pieces for an energy 
eigenstate in Eqs. (9) and (10), along with Eq. (6), imply 



(Bi(E)-S(E)) 3 = 0, 



(12) 



where the boundary conditions were used to eliminate aj 
in favor of aj . The matrix Bf is diagonal in the channels 
and enforces the boundary conditions: 



Bf(E) = (-iy 



z ( cxp(~iK A L) 






exp(— ikbL) 



(13) 



Note that depends on E and K only through the 
relative momenta and Z . 

Equation (12) is of the form of a generalized eigen- 
vector problem. For boxes much larger than the size of 
the interaction, energy eigenstates exist for those energies 
for which the matrix {B^(E) — S(E)*) has an eigenvector 
with eigenvalue zero. Of course, if a matrix has a zero 
eigenvalue, it also has a determinant of zero. Thus, the 
condition for an energy eigenstate for this toy problem 
to fit into a box of width L is equivalent to 



det 



Bl(E)-S(E,Z) =0 



(14) 



Equation (14) should be regarded as a generalized 
Liischcr formula for this multi-channel toy problem. For 
any given 5-matrix one can use it to determine the energy 
levels. Note that the result holds equally well for the case 
where the dynamics are relativistic: the only difference 
is that the Kj are fixed by relativistic kinematics rather 
than non-relativistic kinematics. 

We note that even though Eq. (14) was in a problem 
with only two channels, it immediately generalizes to any 
number of channels. In particular, it also contains the fa- 
miliar one-dimensional two-body single-channel Liischcr 
formula. 



III. A NO-GO THEOREM 

In this section we address the question of whether it 
is possible to fully reconstruct the S'-matrix for the toy 
problem from knowledge of the energy levels in finite 
boxes with periodic boundary conditions. To make the 
problem concrete, we assume that we know the energy 
levels as a function of box size for a finite range of box 
sizes. Given access to only such information about the 
dynamics, can one fully reconstruct the 5-matrix for this 
system at some given energy? As will be shown, the an- 
swer is "no" . 

Suppose that all we knew about the dynamics for the 
two channel problem of the toy model is that it was quan- 
tum mechanical, then the only constraint we have for the 
S-matrix is that it is unitary: S £ U(2) at each energy. 
The most general U(2) matrix can be specified by 4 angle 
variables; a useful parameterization of S is: 



S(E) 



i<t>(E) p i9(E)n(E)-3 



(15) 



where a represents the usual Pauli matrices and the 
energy- dependent unit vector n can be parameterized by 

h = sin(a) cos(/3)i + sin(a) sin(/3)y + cos(a)z, (16) 
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where we have suppressed the energy dependence 
through Eq. (3) of n, a, and f3 for clarity. Thus spec- 
ifying S(E) is equivalent to specifying </>, 9, a and (3 as 
functions of energy; specifying the S-matrix at one en- 
ergy merely requires these functions' values at that en- 
ergy. Note, moreover, that the structure of Eq. (15) is 
such that under the transformation 

(jy^cf) 9^2n-9 (17) 

a — > 7r — a /3— >/3 + 7r 

the S'-matrix is unchanged. This allows us to choose as 
a convention a restricted range of the parameters: < 
(j>< 2ir, < 9 < 2vr, < a < ir and -tt/2 < (3 < tt/2. 

It is easy to show that the energy levels given by the 
generalized Liischcr formula of Eq. (14) are completely 
independent of the parameter (3. This means that it is im- 
possible to determine (3(E) from knowledge of the energy 
levels and hence impossible to determine the S'-matrix. 

To see why the energy levels do not depend on /?, con- 
sider the energy- dependent diagonal matrix 

D(E) = eM*x(E)a 3 ) (18) 

where x(E) is a general energy-dependent function 
that specifies D(E). From the multiplicative prop- 
erty of the determinant and the unitary nature of D, 
it follows that if det [Bf(E) - S(E)] vanishes, then 
dct [Dt(E) (Bl(E) - S{E)) D(E)} does too. Moreover, 
since Bf(E) is diagonal it commutes with D(E). Thus, 
whenever Eq. (14) holds for S-matrix S(E) for some 
given box size L, it also holds for 

S'(E) = D^S(E)D. (19) 

It is easy to see that S'(E) is related to S(E) in the 
following way: 

4>'{E) = 4>{E) 9'(E) = 9(E) (20) 
a'(E) = a(E) P'(E) = /3(E) + 2 X (E) . 

Equation (20) implies that the (3 parameter can be shifted 
by an arbitrary energy-dependent amount without affect- 
ing the energy levels for any box size, provided that the 
boxes are large enough for the generalized Liischer for- 
mula to hold. Clearly the the energy levels alone are 
insufficient to determine (3(E) and hence the S-matrix 
cannot be fully specified. 

It is worth noting that the parameter (3 is associated 
with time-reveral-nonivariant dynamics: for (3 = the S- 
matrix is symmetric but otherwise the S-matrix is gener- 
ically not symmetric. It is well known [28] that the S- 
matrix is symmetric for theories in which the dynamics 
is T-symmetric. 

Thus, while there is a no-go theorem for fully deter- 
mining the S-matrix for the most general case of dy- 
namics, it may appear that it is evaded if one knows 
that the dynamics of the system is T-symmetric since 
for T-symmetric theories one can simply set (3 to zero 



at the outset. However, that is not the case. There 
remains a discrete ambiguity even in the T-symmetric 
case. Suppose that we begin with an S-matrix S(E) 
with f3 = that satisfies the generalized Liischer for- 
mula. Next, we pick D by choosing x(E) to be — tt/2 
for all E. Such a transformation keeps S-matrix sym- 
metric. In that case we know from the argument above 
that S'(E) = D^S(E)D will have identical energies for 
all boxes large enough for the generalized Liischer for- 
mula to apply. Using Eq. (20) and Eq. (17) one sees that 
relationship between the parameters of S'(E) and S(E) 
is given by 

4>{E) = 4>(E) 9'(E)=2-k-6(E) (21) 

a'(E) = 7T - a(E) (3'(E)=0. 

As advertised this is a valid symmetric matrix with f3 = 
which is physically distinct from S(E) provided that 
both 9 and a are nonzero. The nature of this discrete 
ambiguity is that of a sign ambiguity for the off-diagonal 
matrix elements of S. 

To summarize, in the context of this simple toy model 
we have shown that it is not possible to fully specify the 
S-matrix given a complete knowledge of the energy levels 
in finite boxes with periodic boundary conditions. 

While this no-go theorem was derived for this simplest 
of toy models, it is easy to see that holds more gen- 
erally. It trivially generalizes to any 1+1 dimensional 
model with two-body scattering with either rclativistic 
or non-relativistic kinematics and any number of chan- 
nels. In this more general context, the number of dis- 
tinct S-matrices with identical energy levels in a box and 
T- invariant dynamics is 2" -1 where n is the number of 
channels. An analogous no-go theorem in the 3+1 dimen- 
sional case will be more cumbersome to prove because of 
the mixing of partial waves due to box boundary con- 
ditions but there is every reason to believe that such a 
theorem will exist. 

This no-go theorem implies that one cannot fully re- 
construct the S-matrix for multi-channel problems from 
a knowledge of energy levels in finite boxes with peri- 
odic boundary conditions. However, it remains possible 
that interesting and important information about the S- 
matrix may be extracted. For example, in the case of 
T-invariant dynamics, one might hope to fully determine 
the S-matrix with high accuracy modulo the discrete am- 
biguity imposed by the no-go theorem. Such a determi- 
nation would allow one to answer questions such as: given 
incident particles in channel A, what is the probability 
that they will emerge in channel B? In a forthcoming pa- 
per we explore algorithms of this kind and discuss the 
serious technical challenges implementing such an algo- 
rithm for lattice QCD. 
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